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We study the effect of spin-orbit coupling on quantum gates produced by pulsing the exchange 
interaction between two single electron quantum dots. Spin-orbit coupling enters as a small spin 
precession when electrons tunnel between dots. For adiabatic pulses the resulting gate is described 
by a unitary operator acting on the four-dimensional Hilbert space of two qubits. If the precession 
axis is fixed, time-symmetric pulsing constrains the set of possible gates to those which, when 
combined with single qubit rotations, can be used in a simple CNOT construction. Deviations 
from time-symmetric pulsing spoil this construction. The effect of time asymmetry is studied by 
numerically integrating the Schrodinger equation using parameters appropriate for GaAs quantum 
dots. Deviations of the implemented gate from the desired form are shown to be proportional to 
dimensionless measures of both spin-orbit coupling and time asymmetry of the pulse. 



I. INTRODUCTION 

A promising proposal for building a solid-state quan- 
tum computer is based on the notion of using electron 
spins trapped in quantum dots as qubits. 1 In such a de- 
vice, two-qubit quantum gates would be carried out by 
turning on and off the exchange interaction between spins 
on neighboring dots through suitable pulsing of gate volt- 
ages. 

When performing such a quantum gate, if nonadiabatic 
errors 2-4 can be safely ignored, 5 both the initial and fi- 
nal states of the two dots will be in the four-dimensional 
Hilbert space of two qubits. In the absence of spin-orbit 
coupling, and neglecting the dipolar interaction between 
spins, the unitary transformation resulting from such a 
pulsed exchange gate will necessarily have the form 



U = cxp — iXS A ■ Si 



(1) 



where A is a dimensionless measure of the pulse strength. 
This simple isotropic form is a consequence of symme- 
try — if spin and space decouple exactly, as they do 
in the nonrelativistic limit, then the system is perfectly 
isotropic in spin space. Up to an irrelevant overall phase 
the gates (1) are the most general unitary operators with 
this symmetry acting on a two-qubit Hilbert space. 

These isotropic exchange gates are useful for quan- 
tum computation. In conjunction with single qubit rota- 
tions, they can be used in a simple construction of a 
controlled-not (CNOT) gate. 1 It has also been shown 
that, even without single qubit rotations, isotropic ex- 
change gates can be used for universal quantum comput- 
ing with proper encoding of logical qubits. 6 ' 7 

When the effects of spin-orbit coupling are included, 
well-isolated single electron dots will have a two-fold 



Kramers degeneracy and so can still be used as qubits. 
However, when carrying out a quantum gate the total 
spin will no longer be a good quantum number. As a 
result there will inevitably be corrections to the isotropic 
exchange gates (1). Motivated by this fact, a number of 
authors have considered anisotropic gates of the form 

U = exp- 2 A(S A • S B +f3 ■ (S A x S B ) 

+ 1 (S A -S B -0-S A )(P-S B )), (2) 

and shown that they have several useful properties. For 
example, in Ref. 8 it was shown that the CNOT construc- 
tion of Ref. 1 is robust against anisotropic corrections of 
the form appearing in (2). It has also been shown that, 
when combined with a controllable Zeeman splitting, the 
gates (2) form a universal set. 9 

The anisotropic terms which appear in (2) are not the 
most general corrections to (1) which can occur when car- 
rying out an exchange gate in the presence of spin-orbit 
coupling. It is therefore important to ask under what 
conditions these corrections can be restricted to have 
this desired form. The key observation motivating the 
present work is that, up to an irrelevant overall phase, the 
gates (2) are the most general two-qubit quantum gates 
which arc both axially symmetric, i.e. symmetric under 
rotations about an axis parallel to the vector (3 in spin 
space, and symmetric under time reversal (S M — > — S M , 
H = A,B). It follows that if these symmetries can be 
maintained throughout the gate operation, and provided 
nonadiabatic errors can be neglected, the resulting quan- 
tum gate is guaranteed to have the form (2). Of course 
symmetry alone cannot determine the values of A, (3 and 
7. However, in practice we envision these parameters will 
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be determined through experimental calibration rather 
than microscopic calculation. Therefore we emphasize 
symmetry as a useful guiding principle. 

In this paper we study the effect of spin-orbit cou- 
pling on exchange-based quantum gates. For concrete- 
ness we consider a system of two single-electron quantum 
dots in GaAs. The contribution of spin-orbit coupling to 
the exchange interaction between localized spins in GaAs 
has been studied by Kavokin 10 within the Heitlcr-London 
approximation, and by Gor'kov and Krotkov 11 who de- 
rived the exact asymptotic exchange interaction between 
hydrogen-like bound states at large separation. 

Here we follow Ref. 2 and work within the Hund- 
Mulliken approximation, keeping one orbital per dot, and 
allowing double occupancy. In this approximation, the 
effect of spin-orbit coupling is to induce a small spin pre- 
cession whenever an electron tunnels from one dot to an- 
other. The Hamiltonian governing the two-dot system is 
therefore axially symmetric in spin space with the sym- 
metry axis being the precession axis of the spin. If the 
direction of the precession axis does not change while the 
gate is being pulsed, then the resulting quantum gate will 
also be axially symmetric. 

An additional useful symmetry principle, first sug- 
gested in Ref. 12, is that any time-dependent Hamilto- 
nian Hp(t) which is time-reversal symmetric at all times 
t, and which is then pulsed in a time-symmetric way 
(Hp(t) = Hp{—t)) will lead to a gate which can be de- 
scribed in terms of an effective time-independent Hamil- 
tonian H which is also time-reversal symmetric. Here we 
give a proof of this result. 

Taken together these two results imply that, within 
the Hund-Mullikcn approximation, if the spin-orbit pre- 
cession axis is fixed and nonadiabatic errors can be ig- 
nored, the unitary transformation produced by pulsing 
the exchange interaction between two quantum dots will 
necessarily have the desired form (2) provided the gate 
is pulsed in a time-symmetric way. 

This paper is organized as follows. In Sec. II we de- 
rive the Hund-Mulliken Hamiltonian for a double quan- 
tum dot system in the presence of spin-orbit coupling. 
In Sec. Ill we develop an effective spin Hamiltonian de- 
scription which can be applied to pulsing our double dot 
system, and we review the robust CNOT construction of 
Ref. 8. The implications of time-symmetric pulsing are 
then studied in Sec. IV, and in Sec. V we present numer- 
ical results showing the effect of small time asymmetry 
of the pulse. Finally, in Sec. VI we summarize the results 
of the paper. 



II. HUND-MULLIKEN HAMILTONIAN 

We consider a system of two laterally confined quan- 
tum dots with one electron in each dot. For concreteness 
we assume the dots are formed in a two-dimensional elec- 
tron gas (2DEG) realized in a GaAs heterostructure. 



The system is modeled by the Hamiltonian 

H = T + C + H SO - (3) 
Here T + C is the Hamiltonian studied in Ref. 2, where 

h i = ^-(p i --A(r i )) 2 + V(r i ), (4) 
Ira \ c J 

and C — e 2 /e|ri — r2 1 is the Coulomb repulsion between 
electrons. We take the 2DEG the dots are formed in to lie 
in the xy plane, and for GaAs we take m — 0.067m e and 
e = 13.1. For completeness we include a vector potential 
A = (—y,x,0)B/2 which couples the orbital motion of 
the electrons to a uniform magnetic field B = Bz. We 
will see in Sec. Ill that this orbital coupling does not 
affect any of our arguments based on time-reversal sym- 
metry, while a nonzero Zeeman coupling does. 

As in Ref. 2 lateral confinement of the dots is modeled 
by the double-well potential, 

V(x,y) = ^-^(x"-ar+y 2 )- (5) 

This potential describes two quantum dots sitting at the 
points (x,y) — (±a,0). In the limit of large separation 
the dots decouple into two harmonic wells with frequency 

Jo- 
spin-orbit coupling enters the Hamiltonian through 
the term 

H so = J2 n(ki)-Si, (6) 

i=l,2 

where Tik — p — |A. Time-reversal symmetry requires 
that fi(k) is an odd function of k, fi(k) = — fi(— k). 
Thus f2 is nonzero only in the absence of inversion sym- 
metry. 

For definiteness, we take the 2DEG in which the dots 
are formed to lie in the plane perpendicular to the [001] 
structural direction, which then points along the z-axis. 
However, we allow the x-axis, which is parallel to the dis- 
placement vector of the two dots, to have any orientation 
with the respect to the [100] and [010] structural axes. To 
describe the dependence of f2 on k it is then convenient 
to introduce unit vectors ejno] and ejj 10 j which point in 

the [110] and [TlO] structural directions, respectively, and 
define fc[no] = k • e[ 110 ] and fcjjio] = k • e^ 1Q y We then 
have, following Kavokin, 10 

fi(k) ~ (f D - /i?)%no]%io] + (fo + fR)k [Tw] e [lw] . (7) 

Here fo is the Dresselhaus contribution 13,14 due to 
the bulk inversion asymmetry of the zinc-blende crystal 
structure of GaAs, and fp is the Rashba contribution 15 
due to the inversion asymmetry of the quantum well used 
to form the 2DEG. These quantities depend on details of 
the 2DEG confining potential and so will vary from sys- 
tem to system. 
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It was pointed out in Ref. 16 that Hso has a special 
symmetry when fp = ±f R . This can be seen directly 
from (7). When f D = f R (f D = -f R ) the direction of 
fi is independent of k and is fixed to be parallel to e[ 110 ] 
(erj 10 |). The full Hamiltonian (3) is then invariant under 
rotations in spin space about this axis. We will see below 
that this special case has a number of attractive features. 

In the limit of decoupled dots, and ignoring spin-orbit 
coupling, the single electron ground states will be the 
Fock-Darwin ground states centered at (x,y) = (±o, 0), 



Here 



Mx,y) = ^ e -M^a) 2 +y 2 )/2H e ±iay/2ll_ (g) 

V 7ra 

Here u 2 = ui^ + un\ is the frequency of the magnetically 
squeezed oscillator where ujl — eB/2mc is the Larmor 
frequency and l B — yhc/eB is the magnetic length. In 
zero magnetic field, the size of these w ave funct ions is set 
by the effective "Bohr radius" a B = yjTi/mujQ. 

The Fock-Darwin states can be orthogonalized to ob- 
tain the Wannicr states 



y/l - 2Sg - g 2 
1 

^1 - 2Sg - J 2 



(0a ~ g<j>-a), 
(<j)- a - g(j>a), 



(9) 
(10) 



where S = (cf>- a \4> a ) and g = (1 - Vl-S 2 )/S. We 
can then introduce second quantized operators c A (cA a ) 
and c Ba (cB a ) which create (annihilate) electrons in the 
states $^4 and & B with spin a =|, j. 

In the Hund-Mullikcn approximation we keep one or- 
bital per dot and allow for double occupancy. This 
amounts to restricting the full Hilbert space of the prob- 
lem to the six-dimensional Hilbert space spanned by the 
states 



\Si) = -^(c 



t t 



"ai c b]) 



\ S 2) = ^Va[ + c Ll C B T )l°>' 



V2 



|T_) 



L,|o), 



i 



l^o) = -=(c 



t J 



^y-'A^Bl + C Al L LS 



t J 



T )|0>, 



\T+) 



-A1 C B 



T |0). 



(11) 

(12) 
(13) 
(14) 
(15) 
(16) 



In terms of second quantized operators, the Hund- 
Mullikcn Hamiltonian acting in this space, up to an ir- 
relevant overall additive constant, can be written 



Hhm = 



a,/3=T,i 

+V (S A ■ S B + 3/4) 
+U H {n A \n A i +n B \n B [). 



(18) 



<*,/3=U 

is the spin operator on site \x = A, B, 

V = (Si\C\Si) — (T\C\T) (19) 

is the ferromagnetic direct exchange, 

U H = (S 2 \C\S 2 ) - (S^CIS^ (20) 

is the Coulomb energy cost of doubly occupying a dot, 
and 



t H = (® A \h\<f> B ) 



(21) 



is the interdot tunneling amplitude. 

The only contribution from spin-orbit coupling is the 
matrix element 

iP = ($ A |n(k)|* B > = (* A \(p x - e -A x )/n\<z> B ) v , (22) 

where 

V = (Id - fn.) cos 0e [Tlo] + (f D + f R ) sin 8e [lw] . (23) 

Here 9 is the angle the x-axis makes with the [110] struc- 
tural direction. This term introduces a small spin pre- 
cession about an axis parallel to P through an angle 
<f> = 2 arctan(P/i#) when an electron tunnels between 
dots. 

It is convenient to express the spin-orbit matrix ele- 
ment as P = slso where 



v {Id ~ Ir) 2 cos 2 9 + (f D + f R ) 2 sin 2 9 



a B huo 



(24) 



is a dimcnsionless measure of the strength of spin-orbit 
coupling. As stated above, f B and f R depend on de- 
tails of the potential confining the electron to the 2DEG. 
Thus 9, f B and f R are all parameters that in principle 
can be engineered to control the value of s. For example, 
if 9 = then s = \ f B — /ij|/( a B^ w o)- Thus, for this ori- 
entation of the dots, if it is possible to design a system 
in which fn = f R , s can be made to vanish. Even if such 
perfect cancellation cannot be achieved, minimizing the 
difference fo — Ir will reduce s. 

In what follows we leave s as a free parameter. We 
estimate that for GaAs quantum dots s < 0.1 for typical 
parameters. 10 The remaining contribution to the matrix 
element P is then 



(c Aa (t H 8 a p + iP ■ (T a p)c B[j + H.c)j Iso = 



hu! 1 — g 2 d 
~2~ 1 - 2Sg + g 2 b 



e -d 2 K2-l/6 2 )~ ; 



(25) 



where d — a/a B is a dimcnsionless measure of the dis- 
(17) tance between dots, b = yjl + lu l /lo^, and ff = rj/r]. The 
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geometry of our model system is shown schematically in 
Fig. 1. ' 

In what follows we envision pulsing quantum gates by 
varying the distance d between dots as a function of time. 
In doing this, we will assume that throughout the pulse 
the values of fo and /# do not change. If this is the case 
s will be constant and all of the time dependence of P will 
be due to lso- In addition, the direction of the vector P 
will not change as a function of time. The Hamiltonian 
Hrm will therefore be invariant under rotations in spin 
space about a single fixed axis parallel to P throughout 
the pulse. We will refer to such a pulse as having axial 
symmetry. 



B [001],z 




FIG. 1. Sketch of the GaAs double quantum dot system 
considered in this paper. There is one electron per dot, and 
the dot separation is 2a. The dots are taken to lie in the plane 
perpendicular to both the [001] axis and an applied magnetic 
field B. The displacement of the dots makes an angle 8 with 
the [110] axis. Due to spin-orbit coupling electron spins pre- 
cess about an axis parallel to P when tunneling between dots. 

It is important to note that this axial symmetry is ap- 
proximate. In general Jo and fn will depend on time as 
the gate is pulsed, though in principle the system can be 
engineered to minimize this effect. Also, for general fn 
and fn the appearance of only one vector in spin space 
is a consequence of restricting the Hilbert space to one 
orbital per dot. If more orbitals are included then more 
spin-orbit matrix elements will appear in the Hamilto- 
nian, corresponding typically to different spin-precession 
axes, thus breaking the axial symmetry. However, as 
shown above, if fo = then the full Hamiltonian (3) 
is axially symmetric - thus for this special case all spin 
precession axes will be parallel and axial symmetry will 
not just be an artifact of the Hund-Mulliken approxima- 
tion. In Sec. V we discuss the effect deviations from axial 
symmetry will have on our results. 

Given an axially symmetric pulse, it is convenient to 
take the z-axis in spin space to be parallel to P. For this 
choice, the states \T + ) and |T_) decouple, each having 
energy V. 

Another useful symmetry of Hhm is invariance under 
CAa — * CB-a and CBa — * c-A-a- This transformation 
changes the sign of the states \Si), \S-i) and \To), while 
leaving IS3) invariant. It follows that the state IS3) also 
decouples with energy Ur- The matrix representation of 



Hhm in the remaining nontrivial |T ), \Si), \S2) basis is 
then 

/ V -2iP\ 
Hhm = -2t H . (26) 
\2iP -2t H U H J 



III. EFFECTIVE SPIN HAMILTONIAN 

We now consider pulsing the Hamiltonian Hhm by 
varying the distance between the dots, the barrier height, 
or some combination of the two, in such a way that the 
two electron spins interact for a finite period of time, 
but are well separated at the beginning and end of the 
pulse. We assume the initial state of the system is in 
the four-dimensional Hilbert space describing two qubits, 
i.e. the space spanned by the singly occupied states 
|T ), |T_) and \T + ). As the pulse is carried out, the 
eigenstates of Hhm at any given instant in time can be 
grouped into four low-energy states separated by a gap 
of order Uh from two high-energy states. If the pulse 
is sufficiently adiabatic on a time scale set by ~ U/Uh, 
the amplitude for nonadiabatic transitions which would 
leave the system in the excited state IS2) at the end of 
the pulse can be made negligibly small. 5 If this condition 
holds, the final state of the system can also be assumed 
to be in the four-dimensional Hilbert space of two qubits. 
We will see that this condition is easily achieved in Sec. V. 

One way to theoretically study the effect of such a pulse 
would be to first reduce Hhm to an effective anisotropic 
spin Hamiltonian acting on the four-dimensional low- 
energy Hilbert space and then consider pulsing this ef- 
fective model. 12 The problem with this approach is that 
any such effective spin Hamiltonian will only be valid if 
the pulse is adiabatic, not only on the time scale H/Uh, 
but also on the much longer time scale set by the in- 
verse of the small energy splittings within the low-energy 
space due to the spin-orbit induced anisotropic terms. 
However, it is precisely the nonadiabatic transitions in- 
duced by these terms which give rise to the quantum gate 
corrections we would like to compute. 

Although we may not be able to define an instanta- 
neous effective spin Hamiltonian during the pulse, we can 
define one which describes the net effect of a full pulse. 
This definition amounts to parameterizing the quantum 
gate produced by the pulse as 

U = e-* rH , (27) 

where U acts on the four-dimensional Hilbert space of the 
initial and final spin states. H is then an effective spin 
Hamiltonian, i.e. it can be expressed entirely in terms 
of the spin operators and S#, and r is a measure of 
the pulse duration. Note the definition of r is arbitrary 
because it is the product tH which determines U . Here, 
and in the remainder of this paper, we work in units in 
which h = l. 
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If we assume exact axial symmetry throughout the 
pulse, the effective spin Hamiltonian must be invariant 
under rotations about the z-axis in spin space and must 
also leave the states \T + ) and |T_) degenerate. The most 
general such spin Hamiltonian, up to an irrelevant addi- 
tive term proportional to the identity operator, is 

TH(X;a,(3,j) = \(S A ■ S B + ^(S Az - S Bz ) 
+ P(SaxSb v — Sa v Sbx) 

+ -f(SA x SB x +SAyS B y)), (28) 



and we denote the corresponding quantum gate as 
f7(A;a,/3, 7 ) = e -iT»(*<*,Pri. 



(29) 



When a = 0, this is precisely the gate (2) for (3 || z. 

The CNOT construction originally proposed in Ref. 1 
is based on the sequence of gates 



U g = [/(7r/2;0,0,0)e" s ^J7(7r/2; 0,0,0), 



(30) 



where U(ir/2; 0, 0, 0) = exp — i[(n/2)SA • S B ] is a square- 
root of swap gate. The CNOT gate is then 



Ucnot = e ^l^ Aze i{,ms BzU 



(31) 



Remarkably, it was shown in Ref. 8 that if A = tt/2 and 
a = this construction is robust against the (3 and 7 
corrections, i.e. the gate 



U g = U(n/2; 0, (3, 7 )e^ U(n/2; 0, (3, 7) 



(32) 



is independent of (3 and 7. 

For completeness, we briefly review the arguments of 
Ref. 8. Due to axial symmetry, the action of the gate 
U(X;a, (3, 7) on the states \T + ) and |T_) is trivial and 
independent of a, (3 and 7, 



U(\;a,f3, 1 )\T ± )=e-^ 4 \T ± ). 



(33) 



We can then introduce a pseudospin description of the 
remaining space, where \Si) is pseudospin down and |T ) 
is pseudospin up. The action of the gate U (A; a, (3, 7) on 
this pseudospin space is a simple rotation, 



U{\\a,M) => e lA/4 e- 4b - r / 2 , 



(34) 



where b = \{a,(3, 7 + 1) and the components of r = 
(t x ,t V7 t z ) are pseudospin Pauli matrices. At the same 
time, the action of the single qubit rotation entering U g 
is 



(35) 



Thus to show that the CNOT construction is indepen- 
dent of (3 and 7 if a — we need only show that the 
product 



e -ib-T/2 T g-ib-T/2 



(36) 



is independent of (3 and 7 if a = 0. This condition has 
a simple geometric interpretation. It is the requirement 
that a rotation about an axis parallel to b, followed by 
a 180° rotation about the x-axis, and then a repeat of 
the initial rotation must be equivalent to a simple 180° 
rotation about the cc-axis. This will trivially be the case 
if the vector b = A(a, (3, 7 + 1) lies in the yz plane. Thus, 
if a = 0, this condition is satisfied and the CNOT con- 
struction is exact. Conversely, if a ^ the construction 
is spoiled. 



IV. TIME-REVERSAL SYMMETRY 

In this section we prove the following general result. 
Any time-dependent Hamiltonian Hp{t) which is time- 
reversal symmetric for all t, and for which the time de- 
pendence is itself symmetric, i.e. Hp(to — t) = Hp(to+t) 
for all t, will generate a unitary evolution operator U = 
exp —itH where H is a time- independent effective Hamil- 
tonian which is also time-reversal symmetric. We then 
show that this theorem implies that the parameter a, 
which spoils the CNOT construction described in Sec. Ill, 
is equal to zero for time-symmetric pulsing. 

The time-reversal operation for any quantum system 
can be represented by an antiunitary operator O. 17 An 
orthonormal basis {|Mj)} for the Hilbert space of this 
system is then said to be a time-symmetric basis if 



Q\Mi) - I Mi) 



(37) 



for all i. 

For any Hamiltonian H acting on a state |Mj) in this 
basis we can write 



H\M i )=Y,{M j \H\M i )\M j ). 



(38) 



Under time reversal H is transformed into Oi/O" 1 . Us- 
ing the invariance of the {|Mi)} basis and the antiunitar- 
ity of O we can then also write 

BHQ-^Mi) = QH\Mi) (39) 

= QY t {M j \H\M i )\M j ) (40) 
3 

= J2(M j \H\M i )*\M j ). (41) 



Comparing (38) and (41) leads to the conclusion that 
if H is time-reversal symmetric, i.e. H = OiJO -1 , then 
the Hamiltonian matrix is purely real in the {|Mi)} basis, 
while if H is antisymmetric under O, i.e. H = — OiJO -1 , 
then the Hamiltonian matrix is purely imaginary. 

Since H is real in the {|Mi)} basis if and only if H is 
time-reversal symmetric it follows that the unitary oper- 
ator U = exp— irH is self-transpose, i.e. U = U T , if and 
only if H is time-reversal symmetric. 
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Now consider a time-dependent pulse described by the 
Hamiltonian Hp(t). We assume that Hp(t) is time- 
reversal symmetric at all times, i.e. Hp{t) = 6_ffp(t)0 _1 
for all t. The corresponding unitary evolution operator 
U which evolves the system from time ti to tp can be 
written 



U= lim U{t N )U{t N - 1 )---U{t 2 )U{t 1 ) (42) 



where 



U{U) = e 



_ -iAtHp{n) 



(43) 



with At = (tp — tj)/N and t\ = tj and ijv = tp. 

Since Hp(ti) is time-reversal symmetric the above ar- 
guments imply U T (ti) — U(ti) when U(ti) is expressed 
in the time-symmetric basis {|Mi)}. Thus, in this basis, 
we have 



U 1 = lim {U{t N )U{t N - X ) 

N—>oo 



(44) 



= lim U T (t 1 )U T (t 2 )---U T (t N _ 1 )U T (t N ) (45) 

N^OO 

= lim U(t 1 )U(t 2 )---U(t N - 1 )U(t N ). (46) 

N^oo 

For a time-symmetric pulse Hp(ti) = Hp(tM+i-i) and 
so Uifi) = U (fjv+i-i)- This allows us to reverse the order 
of the operators in (46) which then implies 

U T = U. (47) 

Thus if we write U in terms of an effective Hamiltonian, 

U = e- tTH , (48) 

the matrix elements of H must be real in the time- 
symmetric basis. H must therefore be time-reversal sym- 
metric, i.e. H = QHQ- 1 . 

To apply this theorem to the present problem we take 
the time-reversal operator for our two-electron system to 
be 



_ ginSAy gi-xSBy J( 



(49) 



Here the antiunitary operator K is defined so that when 
acting on a given state it takes the complex conjugate 
of the amplitudes of that state when expressed in the 
Hund-Mulliken basis defined in Sec. II. Note that this 
basis is constructed using the Fock-Darwin states, and if 
a magnetic field is present these states will be necessar- 
ily complex valued when expressed in the position basis. 
As defined here, the antiunitary operator K only takes 
the complex conjugates of the amplitudes in the Hund- 
Mulliken basis, it does not take the complex conjugate of 
the Fock-Darwin states themselves. Thus, if a magnetic 
field is present, should be viewed as an effective time- 
reversal symmetry operator. This is a technical point 
which does not affect any of our conclusions (provided 
the Zeeman coupling can be ignored — see below). The 
key property that we will need in what follows is that 



spin changes sign under time reversal, and it is readily 
verified that for our definition of O, 



(50) 



for [i = A, B even in the presence of a magnetic field. 

Under 6, the Hund-Mulliken basis states transform as 
follows, 



e\Si)= \Si) fori =1,2, 3, 

e\T ) = -|T ), 
6|T+)= |T_), 
0|T_)= \T + ). 



(51) 
(52) 
(53) 
(54) 



The states 15^) therefore form a time-symmetric basis for 
the singlet states. A time-symmetric basis for the triplet 
states is given by 



\T ) = i\T 



|T ) = - 7 =(|T + ) + |T_)), 
\f b ) = -j=(|T_) - |T + », 



(55) 
(56) 

(57) 



all of which are eigenstates of & with eigenvalue +1. 

The matrix representation of Hh m in the time-reversal 
invariant |Tq), |Si), IS2) basis is 



V -2P 
Hhm = I -2t H 
-IP -2t H U„ 



(58) 



which is real, reflecting the effective time-reversal symme- 
try of Hhm ■ Note that this would not be the case if Hhm 
included the Zeeman coupling of electron spins to an ex- 
ternal magnetic field. While for typical field strengths 
the Zeeman coupling is small, 2 for some parameters it 
can be comparable to the spin-orbit corrections consid- 
ered here. If this is the case our conclusions following 
from effective time-reversal symmetry will no longer be 
valid. Of course in zero magnetic field exact time-reversal 
symmetry is guaranteed. 

We now consider pulsing a time dependent HnM(t) 
adiabatically so that, according to the arguments of 
Sec. Ill, the resulting gate can be parametrized by an 
effective spin Hamiltonian H. Since at all times t the 
Hund-Mulliken Hamiltonian is time-reversal symmetric, 
if the pulse itself is time symmetric, i.e. HnM{t) — 
HnM(-t) where we take the center of the pulse to be at 
t = 0, then the above theorem implies that the effective 
spin Hamiltonian H will also be time-reversal symmetric. 
Thus H = 6HB- 1 , and since GS^ 1 = -S M this im- 
plies H must be quadratic in the spin operators, and so 
a = 0. The resulting gate will therefore have the desired 
form (2). 

For completeness we also consider here the case of 
time-antisymmetric pulsing. If Hp(t) = —Hp(—t) then 
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FIG. 2. Time dependence of matrix elements appearing in the Hund-Mulliken description of a double quantum dot when 
the displacement of the dots is varied according to (61) with do = 1. Results are for GaAs parameters in zero magnetic field 
with Tiujq = 3 meV and are plotted vs. the dimensionless quantity t/r for two values of the time-asymmetry parameter, r = 
(solid line) and r — 0.1 (dashed line). 



U(t) = e - lAtHp( V = e iAtHp (~V = U(-t)-\ (59) 

and the resulting quantum gate is 

U= lim U(h)U(t 2 )---U(t N/2 ) 

x Uitx/z)- 1 ■ ■ ■ Uih^Uih)- 1 = 1. (60) 

The net effect of any time-antisymmetric pulse is thus 
simply the identity transformation. 

V. MODEL CALCULATIONS 

We have seen from symmetry arguments that time- 
symmetric pulsing of an axially symmetric Hamiltonian, 
such as Hhm when fu and Jr are constant, which is 
itself time-reversal symmetric at all times, will automati- 
cally produce a gate of the form (2), provided the pulse is 
adiabatic so that the initial and final states of the system 
are in the four-dimensional Hilbert space of two qubits. 



It is natural to then ask what the effect of the inevitable 
deviations from time-symmetric pulsing will be on the re- 
sulting gate. To investigate this we have performed some 
simple numerical simulations of coupled quantum dots. 

In our calculations, we imagine pulsing the dots by 
varying the dimensionless distance d between them ac- 
cording to 

d(t)=d 0+ (-±- i y . (61) 

Here do is the distance at the point of closest approach, r 
is a measure of the pulse duration, and r is a dimension- 
less measure of the time asymmetry of the pulse. This 
form describes the generic behavior of any pulse for times 
near the pulse maximum (t = 0). Note that for large \t\, 
and for r ^ 0, the distance d(t) will saturate, and has a 
singularity for negative t. We have taken r to be small 
enough so that the dots decouple long before this leads 
to any difficulty. 
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FIG. 3. Parameters appearing in the effective spin Hamiltonian derived from pulses depicted in Fig. 2. The parameters a, 
13 and 7 are shown as functions of s for the case r = (time-symmetric pulses). For a the quantity a/s is plotted vs. r. We 
have verified that the ratio a/s is essentially independent of s for all values we have considered (|s| < 0.1). 



For our calculations, wc work in zero magnetic field 
and take Tito® = 3 meV and do = 1, corresponding to 
a ~ 20 nm at closest approach. The resulting time de- 
pendences of the parameters in Hum are shown in Fig. 2. 
Note that the spin-orbit matrix element plotted in this 
figure is I so, while the spin-orbit matrix element appear- 
ing in Hhm is P = slsoz where s is the dimensionless 
measure of spin-orbit coupling introduced in Sec. II. 

For a given pulse HuM{t) we integrate the time- 
dependent Schrodinger equation to obtain the evolution 
operator U for the full pulse. If the pulse is adiabatic then 
the matrix elements of U which couple the singly occu- 
pied states |5i) and |T ) to the doubly occupied state 
l^) can be made negligibly small. 5 The quantum gate is 
then obtained by simply truncating U to the 4x4 ma- 
trix acting on the two-qubit Hilbert space. By taking the 
log of this matrix we obtain tH = i log U and thus the 
parameters A, a,/3, 7. Note that when calculating logU, 
there are branch cuts associated with each eigenvalue of 
U, and as a consequence tH is not uniquely determined. 
We resolve this ambiguity by requiring that as the pulse 



height is reduced to zero and U goes continuously to the 
identity that tH — > without crossing any branch cuts. 

We fix the pulse width r by requiring that if we turn 
off spin-orbit coupling (s = 0) we obtain a A = 7r/2 pulse, 
i.e. a square-root of swap. For the parameters used here 
we find this corresponds to taking r = 23.9/wo — 5 ps. 
We have checked that these pulses are well into the adi- 
abatic regime. The magnitudes of the matrix elements 
coupling singly occupied states to the doubly occupied 
state \S 2 ) are on the order of \{S 1 \U\S 2 )\ ~ 1CT 6 and 
\(T \U\S 2 }\~ slO- 6 . 



TABLE I. Symmetry properties of the pulse parameters r 
and s, and gate parameters A, a, [3 and 7 under parity P and 
time reversal T. 
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Once r is fixed, there are two parameters characteriz- 
ing each pulse, s and r, and four parameters character- 
izing the resulting gate, A, a, (3, and 7. The transfor- 
mation properties of these parameters under parity (P) 
and time reversal (T) are summarized in Table I. These 
properties follow from the fact that (i) under time rever- 
sal S M — > — S M and r — > — r, while P = slso^ is invariant, 
and (ii) under parity <-> Sb and P — > P, while r is 
invariant. Note that, as defined in Sec. II, the parameter 
s is positive. Here we allow s to change sign when the 
direction of the vector P is reversed, thus under parity 
s — > — s. 

These symmetry properties imply that if s and r are 
small, the parameters of the effective Hamiltonian will 
be given approximately by 



a ~ C a rs, (62) 

P * Cps, (63) 

7^C 7 s 2 , (64) 

X~X + C x s 2 , (65) 



where the coefficients should be of order 1. For the pulses 
we consider here Ao = tt/2. 

The results of our calculations are shown in Fig. 3. 
Each point corresponds to a separate numerical run. The 
plots for A, (3 and 7 show their dependence on s when 
r = 0. The dependence of the parameter a on pulse 
asymmetry is shown by plotting a/s vs. r. For the s 
values we have studied, up to \s\ = 0.1, the numerical re- 
sults for a/s are essentially independent of s for a given 
r. These results are clearly consistent with the above 
symmetry analysis. 

Now consider carrying out a CNOT gate using the 
scheme reviewed in Sec. III. For this construction to work 
it is necessary that X — tt/2. In our calculations we have 
fixed r so that A = n/2 for s = r = 0. Thus, when 
spin-orbit coupling is included 

A~tt/2 + C a s 2 . (66) 

In order to keep A = tt/2 it will therefore be necessary to 
adjust the pulse width r slightly to correct for spin-orbit 
effects. 

The central result of this paper is summarized by the 
equation 

a ~ C a rs. (67) 

As shown in Sec. Ill, any nonzero a will lead to correc- 
tions to the CNOT construction. For time-symmetric 
pulses r = and these corrections will vanish. Equation 
(67) can then be used to estimate the errors due to any 
time asymmetry of the pulse, and to put design restric- 
tions on the allowed tolerance for such asymmetry. 

It is important to note that while the results presented 
here are for a specific model, all of the key arguments are 
based on symmetry and so are quite general. Given any 



time-reversal invariant two-qubit system with axial sym- 
metry, if pulsed adiabatically in a time-symmetric way 
the resulting gate will have the form (2). 

If the pulse is not axially symmetric, e.g. if the ra- 
tio varies during the pulse, then time-symmetric 
pulsing will still restrict the resulting gate to be invari- 
ant under time reversal. Thus, up to an irrelevant overall 
phase, this gate will necessarily have the form 

U = cxp -iX(S A ■ S B + (3 ■ S A x S B + S A ■ W ■ S B ). (68) 

Here IT is a symmetric tensor which will, in general, de- 
viate from the axial form of the 7 term in (2) leading 
to corrections to the CNOT construction. However, be- 
cause IT is even under parity it will still be second order 
in spin-orbit coupling, 12 and thus the deviations from (2) 
will also be second order. We conclude that even in the 
absence of axial symmetry, the corrections to the CNOT 
construction will be second order in spin-orbit coupling, 
rather than first order. 



VI. CONCLUSIONS 

In this paper we have studied spin-orbit corrections to 
exchange-based quantum gates, emphasizing symmetry 
arguments. In particular, we have shown that adiabatic 
time-symmetric pulsing of any Hamiltonian which (i) de- 
scribes two well defined spin- 1/2 qubits at the beginning 
and end of the pulse, (ii) is time-reversal symmetric at 
all times during the pulse, and (iii) is axially symmetric 
in spin space with a fixed symmetry axis, will automati- 
cally produce a gate of the form (2). Together with single 
qubit rotations, for A = tt/2 this gate can then be used 
in a simple CNOT construction. This result is quite gen- 
eral. 

As a specific example we have studied a GaAs double 
quantum dot system within the Hund-Mulliken approxi- 
mation. In this approximation spin-orbit coupling enters 
as a small spin precession when an electron tunnels be- 
tween dots. If the direction of this precession axis is 
constant throughout the pulse the resulting gate will be 
axially symmetric and have the form (29). The deviation 
of this gate from the desired gate (2) is then character- 
ized by a single dimensionless parameter a which spoils 
the CNOT construction. Using symmetry arguments, 
as well as numerical calculations, we have shown that 
a ~ C a sr where s and r are, respectively, dimensionless 
measures of spin-orbit coupling and time asymmetry of 
the pulse. Thus time-symmetric pulsing (r = 0) ensures 
the anisotropic corrections will have the desired form. 

In any system without spatial inversion symmetry, 
spin-orbit coupling will inevitably lead to anisotropic cor- 
rections to the exchange interaction between spins. Ac- 
cording to current estimates, 18 fault-tolerant quantum 
computation will require realizing quantum gates with 
an accuracy of one part in 10 4 . Thus, even if spin-orbit 
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coupling is weak, the design of any future quantum com- 
puter which uses the exchange interaction will have to 
take these anisotropic corrections into account. We be- 
lieve the symmetry based analysis presented in this paper 
provides a useful framework for studying these effects. 
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